Photoexcited electron dynamics in Kondo insulators and heavy fermions 
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We have studied the photoexcited carrier relaxation dynamics in the Kondo insulator SmB6 and 
the heavy fermion metal YbAgCu4 as a function of temperature and excitation level. The dynamic 
response is found to be both strongly temperature dependent and nonlinear. The data are analyzed 
with a Rothwarf- Taylor bottleneck model, where the dynamics are governed by the presence of 
a narrow gap in the density of states near the Fermi level. The remarkable agreement with the 
model suggests that carrier relaxation in a broad class of heavy electron systems (both metals and 
insulators) is governed by the presence of a (weakly temperature dependent) hybridization gap . 



Recent experiments on the photoexcited electron dy- 
namics in heavy fermion metals have shown that the pi- 
cosecond (ps) relaxation time of the photoinduced (PI) 
reflectivity increases by more than two orders of magni- 
tude upon cooling from 300 to 10 K 0,0] • In contrast, the 
dynamics of their non-magnetic analogues were similar to 
conventional metals like Au and Ag [3j , and could be de- 
scribed by the so-called two-temperature model (TTM) 
0, 0, where the ps recovery is governed by electron- 
phonon (e-ph) thermalization. The relaxation times in 
magnetic and non-magnetic compounds were quite sim- 
ilar at high temperatures (T) suggesting the TTM as a 
starting point for understanding the relaxation processes 
also in heavy fermions. The TTM analysis of the dynam- 
ics showed that the relaxation process in heavy fermions 
can be accounted for using a simple e-ph thermalization, 
assuming that there exists a mechanism for e-ph scatter- 
ing suppression when both the initial and hnal electronic 
states lie within the peak in the density of states (DOS) at 
the Fermi level (E/) 0,0. It was argued that the small 
Fermi velocity compared to the sound velocity may be 
the origin of this suppression, because, in this case, the 
energy and momentum conservation law suppresses e-ph 
thermalization 0- 

In the previous work, there were some observations 
that were not described by the TTM model, such as the 
anomalous rise-time dynamics observed in YbAgCu4 at 
low T 1] . Importantly, similar T-dependent rise-rime dy- 
namics have also been observed in superconductors p|- 
Furthermore, the T-dependence of the relaxation rate ob- 
served in YbAgCu4 also closely resembles the data on su- 
perconductors, both conventional [j| and cuprate [E0- 
The very similar behavior of the PI carrier dynamics 
in superconductors and heavy electron systems suggests 
that the physics governing the relaxation dynamics in 
heavy electron compounds may need to be reconsidered. 
Indeed, if hybridization of the local f-moments with con- 
duction electrons leads to the opening of a well estab- 
lished (indirect) hybridization gap (E g ) near E/ H, 0|, 
similar relaxation bottleneck is expected as observed in 
superconductors . The analogy with superconductors 
is straightforward for the case of Kondo insulators, where 



Eif lies within the gap ||. However, the same bottleneck 
physics 01 is expected to be effective even in the case of 
metallic heavy fermions, if Ef lies close to the hybridiza- 
tion gap edge (if the distance between the gap edge and 

Ef « Eg). 

Relaxation phenomena for nonequilibrium supercon- 
ductors have revealed some of the most intriguing prob- 
lems in condensed matter physics since the 1960's [13 . 
The presence of the superconducting gap in the single 
particle excitation spectrum presents a relaxation bottle- 
neck. The minimal model that describes the relaxation 
of the (photo)excited superconductor was formulated in 
1967 by Rothwarf and Taylor ^lj. Pointing out that the 
phonon channel should be considered when discussing re- 
laxation processes, they described the relaxation dynam- 
ics by two coupled nonlinear-differential equations |ll| . 
While in the low perturbation limit the equations can be 
linearized it was shown only recently that approx- 
imate analytical solutions can be obtained for all limit- 
ing cases [l3j . Analytical solutions enable comparison of 
the experimental data with the model, and have revealed 
that the Rothwarf- Taylor (RT) model can account for 
both the rise-time dynamics [S| as well as the supercon- 
ducting state recovery as a function of excitation fluence 
(F) and T 0. 

In this Letter, we present the first detailed study of the 
carrier relaxation dynamics in heavy electron systems as 
a function of excitation level. Utilizing a low repetition 
rate optical parametric amplifier we were able to mea- 
sure the F -dependence of the reflectivity transients in the 
Kondo insulator SmBg and heavy fermion YbAgCu4 over 
more than three orders of magnitude in F. Both the tran- 
sient amplitude and the relaxation rate were found to be 
strongly F-dependent. Moreover, both observations are 
found to be consistent with the RT model implying that 
the relaxation phenomena in Kondo insulators, as well as 
in heavy fermion metals, are governed by a phonon bot- 
tleneck mechanism arising from the presence of the hy- 
bridization gap in the DOS at Ef. The temperature de- 
pendence of the transient amplitude and relaxation rate 
further supports this scenario, and the extracted values 
of the hybridization gap are in good agreement with the 
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published data. The data suggest that the hybridization 
gap is, to a first approximation, T-independent consistent 
with the periodic Anderson lattice model in the U = 
limit HQ. 

In order to study the independence of the dynamical 
response at low T while minimizing the continuous heat- 
ing of the probed spot due to laser excitation, we used 
an amplified Ti:Al2 03 laser system and an optical para- 
metric amplifier operating at 250 kHz producing sub 100 
fs pulses. The experiments were performed using a stan- 
dard pump-probe set-up. The samples were excited at 
3.0 eV with F ranging from 0.5 to 600 /xJ/cm 2 , while 
the PI changes in reflectivity were measured at a photon 
energy of 1.67 eV. The lowest F used in this study was 
only about a factor of 5 higher than F used in an ear- 
lier study Q, while the maximum F is more than 1000 
times higher. Moreover, since the repetition rate of the 
amplified system is 3 orders of magnitude lower than in 
the high repetition rate system the effect of sample 
heating is minimized. 
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FIG. 1: (a) The F dependence of the PI reflectivity taken in 
YbAgCu4 at 5 K. The data are normalized to F to emphasize 
the sub-linear F-dependence of the amplitude A. (b) The T- 
dependence of the PI reflectivity taken at F — 8.4 /ij/cm 2 . 
The inset to (a) illustrates the way A and r are extracted 
from the raw data; the data set taken at 5 K and F = 0.28 
nJ/cm 2 serves as an example. 



Figure 1(a) presents the PI reflectivity in YbAgCu4 
measured at 5 K while F is varied from 0.28 to 580 
/xJ/cm 2 . All the reflectivity traces in panel (a) are nor- 
malized to F, i.e. AR/R is divided by F in units of 
/ij/cm 2 . For a linear response, all the traces normalized 
to F would lie on top of each other; this is clearly not 
the case. The inset shows the data taken at 0.28 zxJ/cm 2 
on a linear scale. The recovery dynamics is found to be 
non-exponential, and the relaxation non-universal, which 
is to be expected in the non-linear regime. Therefore, 
when performing a quantitative analysis of the data, we 
analyze the amplitude of the transients, A, and the ini- 
tial relaxation rate, r _1 , extracted by fitting the early 



relaxation dynamics with exponential decay (see inset to 
panel (a)). 

At the lowest F, the data are consistent with the low 
T data obtained in an earlier study 0. The initial fem- 
tosecond rise-time (~ 200 fs, resolution limited) is fol- 
lowed by a further increase of the signal on a ps timescale. 
The recovery dynamics proceeded on a 100 ps timescale. 
As F is increased, the PI transient is changed dramat- 
ically. The second-stage risetime becomes faster and 
above F m 4 /iJ/cm 2 the rising edge is resolution lim- 
ited. The decay rate is also weakly independent up to 
w 4 ixJ/cm 2 and increases by almost 2 orders of magni- 
tude as F is further increased to the 100 /xJ/cm 2 range. 
Moreover, the amplitude A shows linear F-dependencc 
up to F « 4 /iJ/cm 2 , while at higher fluences A ps \/~F. 
Panel (b) shows the T-dependence of the reflectivity tran- 
sient taken at constant F. The behavior is similar to 
earlier low F data , revealing two orders of magnitude 
increase in the relaxation rate upon warming to 300 K, 
while A gradually decreases. 

Figure 2 presents F- and T-dependence of the PI re- 
flectivity in the Kondo insulator SmBg |8| . The behavior 
is very similar to YbAgCu4. The main difference is that 
in SmBg no anomalous two-stage risetime is observed. 
The fluence at which A departs from A oc F is near 12 
/iJ/cm 2 . 
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FIG. 2: (a) The independence of the PI reflectivity in the 
Kondo insulator SmBg at 5 K. The data have been normalized 
to F in order to emphasize the sub-linear F-dependence of the 
amplitude, (b) The T-dependence of the PI reflectivity taken 
at F= 8.4 itJ/cm 2 . 



Before analyzing the F- and T-dependence of the PI re- 
flectivity transients, we briefly review the RT model and 
its solutions, using language appropriate for a photoex- 
cited narrow gap semiconductor. Following photoexcita- 
tion, high energy electrons initially release their energy 
via e-e and e-ph scattering. After this process, which 
usually proceeds on a sub-ps timescale 0, the system is 
characterized by excess densities of electron-hole pairs 
(EHP) and high frequency phonons (HFP). When an 
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EHP with an energy > E g (E g is the gap) recombines a 
high frequency phonon (cu > E g ) is created. Since HFP 
can subsequently excite EHP, the recovery is governed by 
the decay of the HFP population . The dynamics of 
EHP and HFP populations is determined by 



dn/dt = nN - Rn 2 , 

dN/dt = -r)N/2 + Rn 2 /2 - j(N - N T ). 



(1) 



Here n and N are the concentrations of EHP and HFPs, 
respectively, n is the probability for EHP creation by 
HFP absorption, and R the rate of electron-hole recom- 
bination with the creation of a HFP. Nt is the concentra- 
tion of HFPs in thermal equilibrium, and 7 their decay 
rate (governed either by anharmonic decay or by diffusion 
out of the excitation volume) . 

The thermal equilibrium concentrations of HFPs and 
EHPs (tit) satisfy the detailed balance equation Rn T — 
(3Nt- Depending on the initial conditions (no and No, 
which are concentrations of EHP and HFP after photoex- 
citation and the initial e-e and e-ph avalanche process) 
and the ratio of 7/77, several different regimes are real- 
ized [13 ■ It follows from the RT analysis that the ps 
risetime dynamics is observed at low F when photoex- 
citation and the initial avalanche processes lead to an 
excess phonon population with respect to the detailed 



balance condition; i.e. Ruq < (3Nq [13|. This is followed 
by the "thermalization" of EHP and HFP distributions 
leading to quasi-stationary distributions of n s and N s 
that satisfy |j| 



Rnj = PN S 



R, IQR 8R„ 

— (4/I+ n Q + — iV 

4n V n n 



(2) 



Since the PI reflectivity (absorption) is proportional to 
the PI carrier density, the resulting PI reflectivity shows 
an initial fast rise-time (n(t) reaches no) followed by a ps 
rise where n(t) reaches n s . Recovery proceeds on a much 
longer timescale and is governed by HFP decay. 

There are several important predictions of the model 
that can be used to determine whether photoexcited car- 
rier dynamics in heavy electron systems indeed follow 
the RT kinetics. First, we should note that the dynam- 
ics at low T should be strongly F -dependent. Since the 
amplitude A is a measure of the PI e-h density, i.e. A 
oc n s — tit, and no and Nq are, at low-T, proportional to 
F, it follows from Eq.@ that A oc y/l + cF — 1, where c 
is a constant. As seen in Fig. 3(a) the independence of 
A for both SmBg and YbAgCu4 is in excellent agreement 
with the model. 

The initial relaxation rate r _1 is given by |l3| 



2R-y(n s 



■ n T 



n 2 (l + 2 7 /n) 



(3) 



It suggests that at low F, when n s k iiy, the relaxation 
rate should be F-independent, while at high fluences, 



r _1 oc F. As shown in Fig. 3(b) the data for SmB 6 
agree well with the model over the entire range of F. In 
YbAgCu4, on the other hand, at high F the experimen- 
tal r" 1 increases faster than the theoretical prediction, 
finally saturating. Saturation is also observed in A(F) 
for both compounds. This saturation can be attributed 
to the PI smearing of the hybridization gap structure. 
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FIG. 3: The independence of (a) the amplitude A and (b) 
the initial relaxation rate r _1 for SmBe and YbAgCu4 at 5 
K. Solid lines are fits to the data with the RT model, A oc 
vT + cF — 1 and r _1 oc (n s + nr) (see text), while dashed 
lines in panel (a) represent the best linear fit (arrows indicate 
the departure from A cc F dependence) . 



The present F-dependence studies clearly support the 
idea of electron relaxation dynamics in heavy electron 
compounds being governed by the presence of a hy- 
bridization gap in the DOS. Moreover, as was shown in 
Ref. 0, the T-dependencies (at constant F) of both A 
and r" 1 are governed by the T-dependence of the num- 
ber density of thermally excited EHP's, jit- It was shown 
that n T oc A^ 1 - 1, where A(T) = A(T)/A(T — ► 0). It 
is easy to demonstrate that the T-dependence of t -1 is 
governed by the T-dependence of nr as well. Assuming 
R, 7, and n to be T- independent, Eq.@ can be rewritten 
as t _1 (T) = C[(n s — ut) + 2ut], where C is a propor- 
tionality constant. Furthermore, since A oc n s — tit and 
tit oc A^ 1 — 1, it follows that 



- 1 (T) = C[D(n T + l)- 1 + 2n T ] 



(4) 



where D is a constant that depends only on the photoex- 
citation intensity. Therefore, for constant F, t _1 (T) is 
entirely determined by the T-dependence of tit- 

In a narrow band semiconductor the T-dependence of 
riT depends on the shape of the DOS in the energy range 
e « kT around the chemical potential. Generally, ut is 
given by: 



n T ~TPexp(-£ g /2T), 



(5) 



where p is on the order of 1, depending on the exact shape 
of the DOS near the gap edge. For a hybridization gap 
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scenario, the shape of the DOS should be close to that 
of a BCS superconductor, where p = 1/2. Neither the 
exact shape of the low energy DOS, nor T-depcndencc 
of Eg is well known in heavy electron systems, therefore 
there will be some ambiguity in determining the precise 
value of the indirect gap E g . However, since the main In- 
dependence in Eq. [5] comes from the exponential term, 
a rough estimate of the size and the T-dependence of E g 
can be obtained. 

In Figure 4(a), we plot the T-dependence of tit for the 
two compounds. In both data sets tit changes by almost 
3 orders of magnitude between 20 and 200 K. Moreover, 
the T-dependence of nx is far from linear (expected in the 
TTM scenario Q) over the entire T-range. Furthermore, 
the absence of any discontinuity or change of slope in tit 
at high-T suggests that the hybridization gap structure 
persists to high T, as observed, for example, in the pho- 
toemission spectroscopy of Kondo insulators 01 • The T- 
dependence of t _1 shown in Fig. 4(b) further supports 
this conclusion: in SmBg, r _1 increases all the way to 
the highest temperatures measured, while in YbAgCu 4 , 
r _1 saturates only above w 160 K. 
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FIG. 4: The T-dependence of (a) the density of thermally 
excited electrons nr (oc A^ 1 — 1), where A(T) is plotted in 
the inset, and (b) the initial relaxation rate, r _1 . The data 
were taken at F = 8.4 /zJ/cm 2 . The dashed lines are fits to 
the data by the RT model (see text). 

Based on these qualitative findings, we can apply the 
RT model to extract the magnitude of the indirect hy- 
bridization gap, Eg. We consider the hybridization gap 
to be present at all T, and to a first approximation, to 
be T- independent. Furthermore, we assume a BCS-like 
DOS near the gap edge, i.e. p = 0.5. Fitting tit with 
Eq. ||SJ,we find, for SmBg, an excellent agreement over 
the entire T-range with E g ss 550 K. In YbAgCu4, on 
the other hand, the agreement is good up to « 120 K, 
with the extracted value of E g » 100 K. At higher T in 
YbAgCu4, tit starts to increase faster possibly due to a 
partial suppression of the gap at high T. 

Qualitative agreement between the relaxation rate 
data and the model described by Eq.© is also straight- 



forward. At the intermediate temperatures, n s — hit <C 
tit, and r _1 is governed by the T-dependence of tit ex- 
hibiting cxp(— E g /2T) behavior. At low enough T, how- 
ever, n s 3> tit and the relaxation time saturates, as ob- 
served experimentally. We fit t -1 (T) data with Eq. Q 
and find remarkable agreement. The extracted values of 
the gap are E g sa 85 K for YbAgCu 4 and E g « 350 
K for SmBg, somewhat lower than the values extracted 
from the fit to tit. This can be partially attributed to 
the fact that the T-dependence of R, 7, and rj was ne- 
glected, as well as to the microscopic details like gap 
anisotropy. Comparison with literature shows that for 
SmBg there is a large spread of published values for E g 
|l5|: however, our estimate is in close agreement with 
the T-independent pseudogap energy scale of 290 cm -1 
(450 K) from a recent Raman scattering study 15] . The 
value of Eg 100 K for YbAgCyuis also in close agree- 
ment with recent optical data |l6|]. The overall agree- 
ment and self-consistency of the data with the RT model 
presents a strong argument that relaxation kinetics in 
heavy electron systems is indeed governed by the pres- 
ence of a weakly T-dependent hybridization gap, as for 
a peri odic Anderson lattice model in the U — limit 

We have presented studies of photoexcited electron re- 
laxation and recombination dynamics in the Kondo in- 
sulator SmBg and heavy fermion YbAgCu 4 as a func- 
tion of temperature and excitation level. The dynami- 
cal response is described well using the phenomenologi- 
cal RT model, suggesting that the carrier relaxation and 
recombination dynamics in heavy electron systems are 
governed by the presence of a weakly T-dependent hy- 
bridization gap. 
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